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“A 、 9 -4 一 一 
流体 计算 月 于 


流体 计算 是 当代 计算 数学 、 流 体力 学 和 计算 机 科学 相 结合 的 产物 ， 是 一 门 具 有 强大 生命 力 
的 边缘 科学 。 它 以 电子 计算 机 为 工具 ， 应 用 各 种 离散 化 的 数学 方法 ， 对 流体 力学 领域 的 各 类 问 
题 进行 数 值 实验 、 计 算 机 模拟 和 分 析 研 究 ， 以 期 解决 各 种 复杂 的 工程 实际 问题 。 为 推动 计算 数 
学 与 应 用 数学 理论 、 方 法 的 研究 ， 活 跃 学 术 气 氛 ， 繁 荣 我 国 流 体 计算 的 学 术 研 究 ， 最 大 限度 地 
发 挥 数 学 理论 方法 在 解决 流体 计算 以 及 其 它 学 科 、 各 种 工程 与 实际 应 用 领域 中 的 问题 的 作用 ， 
本 刊 特 组 建 该 期 有 关 流 体 计 算 的 专辑 。 

本 专辑 共 收 录 来 稿 8 篇 ， 可 分 为 3 种 类 型 : 类 型 I 是 关于 有 较 强 应 用 背景 的 数值 算法 研究 
的 ， 有 3 篇 文章 ， 详 见 前 3 篇 文章 ， 类 型 开 主要 是 关于 数值 模拟 计算 的 ， 有 3 篇 文章 ， 详 见 
第 4 至 第 6 篇 文章 ;类 型 II 是 关于 数值 分 析 和 数值 计算 的 ， 有 2 篇 文章 ， 详 见 第 7 和 第 8 篇 文 
佛 。 

在 类 型 1 的 第 1 篇 文章 中 ， 王 晓 东 等 考虑 了 采用 基于 Doi 理 论 的 微 -宏观 双 尺 度 模型 ， 研 究 了 
液晶 聚合 物 在 均匀 剪 切 流 场 中 的 微观 结构 。 为 保证 模拟 结果 的 可 靠 性 ， 在 模型 求解 中 率先 使 用 
了 高 精度 、 高 稳定 性 的 无 网 格 方法 ， 并 通过 对 均匀 流 场 中 均 质 液晶 聚合 物 系 统 的 数值 模拟 ， 
验证 了 该 方法 的 有 效 性 。 文 中 着 重 研究 了 De 数 对 平板 Couette 流 中 液晶 聚合 物 微观 结构 的 影 
响 。 模 拟 得 到 了 液晶 聚合 物 分 子 四 种 典型 的 微观 运动 结构 ， 即 弹性 稳 态 、 周 期 翻转 、 周 期 摆动 
和 随 流 取向 。 定 量 分 析 了 De 数 对 两 种 周期 运动 的 影响 ， 预 测 了 各 种 结构 下 分 子 的 有 序 程 度 和 
可 能 出 现 的 缺陷 ， 并 通过 对 取向 角 的 详细 分 析 ， 发 现 了 两 种 边界 层 效 应 。 

在 类 型 I 的 第 2 篇 文章 中 ， 李 娴 娟 、 许 传 炬 考虑 了 描述 在 离子 浓度 梯度 VC 及 电场 VV 共 
同 存在 的 情况 下 ， 穿 过 渗透 膜 的 离子 (如 钙 ， 钾 ， 钠 ， 氯 ， 镁 等 离子 ) 流 了 了 的 Nernst-Planck 方 
程 。 然 而 ， 计 算 Nernst-Planck 方 程 的 数值 解 会 遇 到 一 些 困 难 。 首 先 ， 由 于 Nernst-Planck 方 
程 是 多 维 的， 那么 对 其 数值 求解 比较 费时 。 其 次 ， 在 实际 问题 中 边界 值 是 未 知 的 ， 要 确定 
合理 的 边界 条 件 并 不 简单 。 目 前 已 有 一 些 关 于 整数 阶 Nernst-Planck 方 程 数值 计算 的 研究 工 
作 ， 提 出 了 时 间 方 向 用 一 阶 差 分 、 空 间 方向 用 中 心 差分 对 Nernst-Planck 方程 进行 离散 ， 用 
有 限 差分 /有 限 元 法 求解 扩展 的 Poisson-Nernst-Planck 方 程 。 但 是 分 数 阶 方程 的 求解 与 整数 阶 
相 比 具有 本 质 的 不 同 。 该 文 考虑 用 以 描述 神经 细胞 中 离子 反常 扩散 现象 的 电缆 型 简化 的 分 数 
阶 Nernst-Planck 方 程 ， 提 出 一 个 时 间 有 限 差分 /空间 谱 元 法 对 该 方程 进行 数值 求解 。 作 者 给 出 
了 数值 方法 的 详细 构造 过 程 以 及 实现 方法 。 数 值 结果 表明 数值 解 在 空间 方向 上 具有 指数 阶 收敛 
精度 ， 在 时 间 方 向 上 具有 2 一 w 阶 精度 。 最 后 ， 通 过 计算 一 个 具有 实际 背景 参数 的 问题 说 明 所 
提 方 法 的 潜在 应 用 。 

新 型 飞行 器 的 气动 设计 和 研制 要 求 对 具有 分 离 、 旋 涡 、 激 波 / 附 面 层 干 扰 等 复杂 流动 特性 的 
气动 问题 采用 雷诺 平均 Navier-Stokes (RANS) 方程 ， 求 解 这 类 方程 往往 要 求 计 算 达 到 千 万 网 格 
点 的 规模 。 对 于 这 种 大 规模 高 精度 的 流体 力学 计算 问题 ， 人 们 常 利 用 高 性 能 计算 机 进行 并 行 计 
算 以 提高 计算 速度 。 并 行 计算 作为 CFD 研究 的 一 个 主要 方向 ， 得 到 了 国内 外 研究 人 员 的 广泛 
重视 和 快速 发 展 。 在 类 型 I 的 第 3 篇 文章 中 ， 郑 秋 亚 等 讨论 了 基于 连续 拼接 多 块 结构 化 网 格 ， 
通过 求解 雷诺 平均 Navier-Stokes 方 程 研 究 并 行 计算 中 的 负载 平衡 问题 。 利 用 组 合 优化 中 的 排 
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序 理 论 设计 负载 下 衡 算法 ， 实 现 了 网 格 数 据 的 自动 划分 和 各 处 理 机 上 计算 任务 的 自动 分 配 。 
在 工作 站 集群 MPI 并 行 环境 下 ， 通 过 实例 考察 了 负载 平衡 算法 和 并 行 计算 的 性 能 ，16 处 理 机 
上 的 负载 均 方 美和 负载 相对 均 方 状 分 别 为 0.0084 和 0.1347 台 ， 并 行 计 算 结果 和 实验 数据 吻合 良 
好 ， 并 行 效率 高 。 该 文 算法 具有 良好 的 可 扩展 性 ， 适 用 于 MIMD 结构 计算 机 上 基于 多 块 结构 
化 网 格 并 行 计算 中 的 负载 平衡 问题 。 

在 类 型 开 的 第 1 篇 文章 中 ， 王 坪 、 田 振 大 建立 了 一 种 基于 投影 法 的 求解 不 可 压缩 Navier- 
Stokes (N-S) 方程 的 高 精度 紧 致 益 分 格式 。 该 方法 时 间 上 采用 Kim 和 Moin 二 阶 投影 法 离散 、 
空间 上 采用 高 精度 紧 致 格式 离散 ， 并 提出 了 一 种 新 的 离散 压力 边界 的 紧 致 格式 ， 同 时 对 
计算 结果 进行 分 析 以 验证 该 投影 法 精度 和 格式 稳定 性 。 文 中 的 Taylor 涡 列 数值 计算 结果 表 
明 ，Kim 和 Moin 投 影 法 能 使 得 压力 场 和 速度 场 均 达到 时 间 : 二 阶 精度 ， 且 高 精度 紧 致 格式 投影 
法 也 具有 空间 高 阶 精度 。 了 驱动 方 腔 数 值 模拟 结果 显示 ， 该 文 对 N-S 方程 的 离散 格式 具有 很 好 的 
可 靠 性 ， 适 用 于 对 复杂 流体 流动 的 小 尺度 问题 的 数值 模拟 和 研究 。 

在 类 型 本 的 第 2 篇 文章 中 ， 尚 月 强 、 何 银 年 基于 完全 重 受 型 区 域 分 解 技 巧 ， 提 出 了 一 种 求 
解 非 定 常 Stokes 方程 的 有 限 元 并 行 算法 。 该 算法 的 基本 思想 是 首先 对 空间 施行 完全 重奏 型 区 域 
分 解 ， 然 后 各 个 处 理 器 使 用 向 后 Euler 格式 独立 并 行 求解 关于 时 间 的 常 微分 方程 ， 在 整个 关于 
时 间 的 迭代 过 程 中 ， 无 需 处 理 器 间 的 通信 ， 具 有 良好 的 并 行 性 能 。 该 算法 中 每 个 处 理 串 所 负责 
的 子 问 题 是 一 个 全 局 问题 ， 它 定义 在 整个 求解 区 域 上 ， 但 绝 大 部 分 白 由 度 米 自 其 所 负责 的 子 区 
域 ， 从 而 使 得 该 算法 稍 加 修改 现 有 的 串 行程 序 即 可 实现 相应 的 并 行 计算 ， 实 现 简单 ， 具 有 重要 
的 使 用 价值 。 同 时 通过 数值 算 例 ， 在 曙光 集群 并 行 机 上 编程 实现 了 上 述 算法 ， 验 证 了 其 有 效 
性 

在 类 型 工 的 第 3 篇 文章 中 ， 段 献 葆 、 秦 新 强 将 经 典 的 形状 关 敏 度 分 析 方法 与 一 种 改进 
的 水 平 集 方法 相 结合 ， 给 出 了 Navier-Stokes 问 题 形 状 优化 的 一 种 新 方法 。 该 算法 是 在 固定 
的 Buler 网 格 上 进行 计算 且 在 优化 过 程 中 不 需要 对 水 平 集 函数 进行 重新 初始 化 ， 从 而 可 以 有 效 
地 节省 计算 时 间 。 数 值 算 例 说 明 该 算法 是 稳定 、 高 效 的 。 

在 类 型 II 的 第 1 篇 文章 中 ， 王 爱 文 等 针对 低 阶 协调 有 限 元 对 @1 - 瓦 , RR -~ PD， 对 二 维 定 
常 不 可 压缩 Navier-Stokes 方 程 ， 提 出 了 建立 在 局 部 压力 投影 上 的 一 类 稳定 化 有 限 元 方法 。 对 
于 稳定 化 有 限 元 逼近 解 (wn,pn)， 得 到 了 该 方法 的 最 优 误差 估计 。 并 做 了 若干 数值 实验 ， 表 明 
了 该 方法 有 很 好 的 稳定 性 和 收敛 性 。 

在 类 型 III 的 第 2 篇 文章 中 ， 穆 保 英 等 基于 变 分 多 尺度 方法 考虑 了 非 定常 Navier-Stokes 方 各 
的 一 种 稳定 化 方法 ， 并 利用 一 个 与 实际 误差 等 价 的 后 验 误差 估计 子 ， 结 合 自 适应 算法 ， 得 到 了 
非 定常 Navier-Stokes 方 程 的 自 适应 变 分 多 尺度 稳定 化 方法 。 证 明了 该 方法 的 稳定 性 及 后 验 整 
体 误差 与 先 验 误差 的 等 价 性 。 该 方法 简单 、 直 观 ， 且 易于 程序 实现 。 最 后 作者 给 出 了 经 典 的 方 
腔 流 数值 算 例 ， 验 证 了 该 方法 的 有 效 性 。 

描述 流体 运动 规律 的 偏 微分 方程 是 一 种 典型 的 非 线性 发 展 方程 ， 由 于 人 们 对 非 线性 现象 本 
质 的 认识 有 限 ， 因 而 数值 模拟 就 成 为 一 种 十 分 重要 的 研究 手段 。 但 直接 数值 模拟 流体 运动 一 个 
很 大 的 困难 就 是 巨大 的 解 题 规模 、 长 时 间 积 分 运算 与 有 限 计算 资源 及 算法 稳定 性 之 间 的 矛盾 ， 
因此 ， 构 造 和 研究 具有 良好 性 能 的 流体 计算 方法 就 显得 尤为 重要 。 这 次 出 版 流体 计算 专辑 是 我 
们 推动 流体 计算 发 展 的 初步 尝试 ， 今 后 我 们 将 会 继续 关注 国内 流体 计算 的 发 展 ， 收 集 和 出 版 更 
新 的 流体 计算 研究 成 果 ， 以 便 引 起 更 多 的 有 关 科研 人 员 关注 流体 计算 的 建设 、 发 展 和 提高 。 


何 银 年 
2010 年 2 月 于 西安 交通 大 学 
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附 : 特 邀 编辑 简介 

何 银 年 : 男 ，1953 年 7 月 出 生 ， 西 安 交 通 大 学 三 级 教授 ， 博 士 生 导师 ， 享 受 政府 特 殊 津 
贴 专家 ， 新 疆 大 学 “天 山 学 者 ”讲座 教授 。1982 年 1 月 毕业 于 陕西 师范 大 学 数学 系 ， 获 学 士 
学 位 ; 1986 人 年 3 月 毕业 于 西安 交通 大 学 数学 系 ， 获 理学 硕士 学 位 : 1992 年 12 月 毕业 于 西安 
交通 大 学 能 源 系 ， 获 工学 博士 学 位 。 博 士 毕业 后 一 直 在 西安 交通 大 学 理学 院 从 事 教 学 和 科 
研 工 作 。 现 担任 陕西 省 计算 数学 学 会 理事 长 、 陕 西 省 数学 学 会 常务 理事 、 中 国 计 算 数学 学 
会 常务 理事 ， 是 两 种 国际 期 刊 《Advances in Numerical Analysis》 和 《Advances in Applied 
Mathematics and Mechanics》, 以 及 三 种 国内 杂志 《数值 计算 与 计算 机 应 用 》,《 计 算 物理 》 和 

《高 校 计算 数学 学 报 》 的 编 委 。 曾 应 邀 赴 和 荷兰 、 美 国 、 加 拿 大 、 香 港 等 地 多 所 大 学 进行 学 术 
交流 。2004 人 年 6 月 在 第 四 届 中 - 瑞 计 算数 学 国际 会 议 作 邀 请 报告 ， 2008 年 8 月 在 中 -日 - 韩 三 边 
计算 数学 国际 会 议 作 邀请 报告 。 已 主持 国家 自然 科学 基金 项 目 4 项 ， 作 为 骨干 和 课题 副 组 长 
分 别 参 加 了 973、863 计 划 项 目 各 1 项 。 研 究 成 果 获 国家 自然 科学 二 等 奖 1 项 ， 省 部 级 自然 科 
学 和 科学 技术 一 、 二 等 奖 共 3 项 。 何 教授 现 已 在 国内 外 顶级 杂志 ， 如 《Numer Math》,《Math 
Comp》,《Comput， Methods Appl. Mech. Engrg.》,《SIAM J Numer Anal》,《IMA J. Numer. 
Anal》 等 发 表 学 术 论 文 123 篇 ， 其 中 被 SCI 收 录 87 篇 ， 被 SCI 文 章 他 引 92 篇 次 。 

近年 来 ， 何 教授 在 与 流体 计算 相关 的 研究 领域 所 取得 的 主要 学 术 成 就 包括 : 

1) 提出 了 求解 粘性 不 可 压缩 流体 运动 问题 的 时 - 空 两 层 或 多 层 算法 (包括 有 限 元 方法 和 
谱 方法 )， 采 用 时 - 空 两 层 或 多 层 算法 求解 非 定常 粘性 不 可 压缩 流体 运动 时 ， 在 时- 空 粗 网 格 上 
求解 非 线性 Navier-Stokes (N-S) 方程 ， 在 时 - 空 细 网 格 上 求解 线性 化 N-S 方程 (Stokes 方程 )， 当 
时 - 空 粗细 网 格 尺度 保持 某 种 匹配 关系 时 ， 所 得 结果 与 在 时 - 空 细 网 格 上 求解 非 线 性 N-S 方程 时 
具有 同样 的 收敛 速度 和 相同 的 稳定 性 ， 这 样 一 来 对 于 求解 长 时 间 发 展 的 N-S 方程， 具有 节约 计 
算 量 的 优点 ， 其 中 空间 离散 可 以 用 有 限 元 方法 或 谱 方 法 。 此 外 ， 我 们 将 两 层 有 限 元 方法 成 功 应 
用 到 了 求解 定 第 的 粘性 不 可 压缩 流体 运动 。 

2) 将 求解 定常 线性 Stokes 问 题 的 稳定 化 有 限 元 方法 应 用 到 研究 定常 和 非 定 常 非 线 
性 Navier-Stokes (N-S) 方程， 探讨 了 方法 的 稳定 性 ， 收 敛 性 和 数值 计算 等 一 系列 问题 。 

3) 研究 了 求解 非 定常 N-S 方 程 的 全 离散 加 宰 有 限 元 方法 的 稳定 性 和 收敛 性 ， 得 到 了 最 优 
的 依赖 于 时 间 步 长 和 网 格 尺度 的 误差 估计 。 

4) 以 往 的 数值 分 析 研究 表明 求解 非 定 常 N-S 方 程 的 一 阶 和 高 阶 隐 / 显 式 格式 ， 要 求 有 稳定 
性 条 件 才 能 保证 解 的 稳定 性 和 收敛 性 ， 即 要 求 时 间 离 散步 长 选取 依赖 于 网 格 离散 尺度 h， 这 
样 使 得 在 每 个 时 间 步 计算 时 ， 只 能 取 小 的 时 间 步 长 。 通 过 重新 分 析 现 有 格式 的 稳定 性 和 收敛 
性 ， 得 出 了 对 于 光滑 的 初始 数据 ， 求 解 非 定常 N-S 方 程 的 隐 / 显 式 格 式 是 几乎 无 条 件 稳定 和 收 
敛 的， 从 而 保证 了 求解 非 定常 问题 可 以 采取 大 时 间 步 长 ， 加 快 计算 速度 ， 具 有 节约 计算 量 的 优 
点 ， 这 对 于 研究 N-S 方程 解 的 长 时 间 行 为 有 非常 重要 的 意义 。 
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